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Abstract. 

We investigate the counting statistics of non-thermal equilibrium, or rare, quantum 
jump trajectories associated to Markovian exciton transport. By applying the large 
deviation theory to open quantum system dynamics, we are able to uncover dynamical 
transitions, or crossovers, in the statistics of such rare trajectories. We find that the 
crossovers are activated by electronic coherences, and occur in systems with a minimum 
of three interacting molecules. The formalism thus reveals subtle consequences of non- 
equilibrium thermal fluctuations on the steady-state statistics of systems exhibiting 
Markovian transport. 

1. Introduction 



Full counting statistics (FCS) is a powerful mathematical tool that enables the 
characterization of a counting process by systematic computation of the generator 
of its statistics. The technique was originally developed for theoretical modelling 
of mesoscopic charge transport [U E], where current fluctuations can provide vital 
information on the nature of the transport mechanism [31 IU [51 El [7J |8]. FCS of a 
classical, or an open quantum system, at thermal equilibrium can be obtained by 
identifying the generating function of the statistics from the corresponding master 
equations [911101 HI]. However, due to the environmental perturbations, non-equilibrium 
fluctuations can easily arise in open quantum systems. This raises the question that 
how are the counting statistics altered by such perturbations, and whether a systematic 
theory of non-equilibrium FCS for open quantum systems can be formulated. 
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The difficulty in addressing this problem lies in the different nature of the two in- 
fluences: thermal fluctuations of the environment are of a classical character and are 
described via the tools of non-equilibrium thermodynamics. The interaction of the quan- 
tum system with the environment is described via the theory of open quantum systems 
[12] . There has been a historical distinction between the two approaches, rendering the 
study of their mutual influence a non-trivial one. 

An important connection between the two descriptions was provided by Garrahan 
and Lesanovsky in a recent contribution [13]. The link between the two descriptions 
was provided through the well-known mathematical theory of large deviation principle 
[13]. In statistical mechanics, large deviation principle provides a framework to deduce 
thermodynamic quantities from sub-ensembles of particles in the configuration space. 
In the theory of open quantum systems, however, it serves as a tool for statistical anal- 
ysis of rare (i.e.: non-equilibrium) quantum-jump trajectories associated to Markovian 
dynamics. A quantum-jump trajectory is a realization of the system's dynamics sub- 
jected to particular bath-induced measurements at stochastic intervals, and represents 
a sequence of configurations in the time domain. By incorporating a coupling field, or 
a dynamical order parameter, that quantifies perturbations from equilibrium, the large 
deviation approach grants access to the complete statistics of the trajectories, both in 
and out of equilibrium. At equilibrium, large deviation theory recovers the results of 
the full counting statistics as expected. Garrahan and Lesanovsky coined the term the 
s-ensemble theory to a technique of mapping out the out of equilibrium trajectories of a 
quantum master equation via the large deviation approach. This terminology will also 
be used here. 

In the context of open quantum systems, the s-ensemble approach allows the iden- 
tification of statistical transitions in the time domain, thereby providing novel physical 
insight on the dynamics of quantum systems subjected to Markovian noise. Dynamical 
transitions have so far been predicted in photon counting [13] [15l [16], [T7] , electron count- 
ing [18], and quadrature measurements [19], suggesting the possibility of their existence 
in a host of other systems. The significance of transitions in the statistics of quantum 
processes, or the means by which they may be observed, are not always clear and tend to 
be problem-specific. It is evident, however, that such novel dynamical features cannot 
be identified from the analysis of the density matrix operator alone. 

In this article we apply the s-ensemble theory to investigate the statistics of rare 
trajectories of quantum evolutions describing excitation energy transport in a molecular 
complex. Assuming strong electronic interactions, and weak system-environment cou- 
pling, exciton relaxation can be described by Lindblad type evolution [12] . The exciton 
relaxation process results in incoherent transfer between electronic eigenstates, reaching 
a steady-state where the average population of the eigenstates obeys the thermal Boltz- 
mann distribution. At a microscopic level, however, transfer of excitation still occurs in 
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a dynamical manner in the steady-state. For instance for a system of two eigenstates at 
thermal equilibrium, detailed balance dictates that ki^Pi — ^2->i-P2 where k^j is the 
rate of energy transfer from the state % to the state j and Pi is the thermally-equilibrated 
population of the state i. Out of equilibrium fluctuations result in deviations from the 
detailed balance. Therefore, in the present context, the s-ensemble theory serves as a 
mathematical tool to investigate changes in the counting statistics as the steady-state 
deviates from thermal equilibrium. We demonstrate that despite the classical nature 
of the transitions, the statistics of the trajectories are not, in general, Poissonian. We 
observe temperature-dependent dynamical transitions for systems containing a mini- 
mum of three sites, and properties such as bunching and anti-bunching that are often 
associated with non-classical correlations. 

This article is organized as follows: in section |2] we introduce the conceptual and 
mathematical framework of the large deviation principle and discuss its application 
to counting statistics. In section [3] the Markovian Lindblad master equation and the 
s-ensemble theory are introduced. In section H] we apply the s-ensemble theory to 
investigate the statistics of rare trajectories in systems of two, three, four and seven 
sites and present numerical results that demonstrate the dynamical transitions in the 
statistics of trajectories. The electronic and vibrational parameters of the Fenna- 
Matthews-Olson (FMO) photosynthetic complex are used in the computation of the 
numerical results, as the Lindblad master equation considered in this article has been 
found to make a reasonable prediction of the steady-state of the complex. Section [5] 
summaries the main conclusions of the article. 

2. Large deviation principle and counting statistics 

A counting process can be decomposed into a number of trajectories where in each 
trajectory the process occurs at different stochastic intervals [20]. In the context of 
quantum jumps, the counting process is the number of jumps observed after a given 
time t. Let us denote the probability of observing K jumps after a time t by P t (K), 
where Y^k=o^(-^) = 1' an( ^ — Pt{K) < 1}- F° r a counting process that satisfies the 
large deviation theory, Pt{K) acquires a large deviation form for large t. That is 



where the function <p(k) is known as the rate function [T3], and k = K/t is the jump 
rate. For a simple (unimodal) counting process, tp(k) has a local minima corresponding 
to the most probable jump rate after a long observation time. The generating function 
of a given probability distribution is define to be 



Pt(K) ~ e 



ttp(k) 



(1) 




(2) 



K=0 
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where s is a conjugate variable to K and is known as the activation parameter. If the 
counting process has a large deviation form, the generating function can be written as 

Z t (s) ~ e"«. (3) 

Z t (s) is now referred to as the large deviation function, while 9(s) is the cumulant 
generating function. The cumulant generating function and the rate function are related 
via a Legendre transform, 

9(s) = - mm[ip(k) + ks] (4) 

k 

The average and the standard deviation of the number of jumps at equilibrium can be 
deduced from the cumulant generating function, and are given by 

(k) = -9'(0) (5) 

(fc 2 > - (fc> 2 = 0"(O). (6) 

Similarly, higher derivatives of 9(s) evaluated at s — 0, are the higher order statistical 
cumulants. Full counting statistics (or the equilibrium statistics) is therefore recovered 
by looking at the derivatives of the cumulant generating function at s = 0. The 
applicability of the large deviation approach to non-equilibrium processes becomes 
evident when the behaviour of the cumulant generating function away from s = 
is considered. To illustrate this let us construct an associated stochastic process defined 
by the probabilities [21] 

q t (K,s) = -± r -P t (K)e- sK (7) 
Z t (s) 

For s = we recover the original counting process. For s < (s > 0) the exponential 
factor amplifies (lowers) the probabilities of rare events. The coupling field s, thus 
provides a systematic means by which rare events can be accessed: unlikely events of 
the counting process {P t (K)}, correspond to typical events of the associated process 
{q t (K,s)}. Rare events can thus be accessed without generating a large number of 
typical events and relying on the statistics of the original counting process. In non- 
equilibrium thermodynamics, the coupling field s is interpreted as a means by which the 
system is perturbed out of equilibrium. The influence of non- equilibrium perturbations 
are therefore deviations from the thermodynamic probabilities. A rare event corresponds 
to a particular trajectory of the system that has been significantly perturbed out of 
equilibrium. 

In statistical mechanics the counting process is the number of microstate 
configurations accessible to the system under certain constraints, making Z t (s) the 
partition function. Analogously, for an open quantum system, Z t (s) is a measure of 
the number of trajectories that can participate in the dynamics of a particular non- 
equilibrium configuration. One can also identify tp(k) as the entropy, and 9(s) as the 
free energy. Extrema of (p(k) therefore correspond to the equilibrium configurations, or 
the typical trajectories. The variable s is (minus) the chemical potential s = — /i [21J. 
Free energy serves as an order parameter: singularities in free energy predict first-order 
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thermodynamic phase transitions, or dynamical transitions in the space of trajectories 
of a quantum system. 

3. s-ensemble theory and electronic energy transfer 

In this section we apply the large deviation theory to investigate the trajectories of a 
particular open quantum system, modelled by the Markovian Lindblad master equation. 
The system of interest consists of a number of molecules (sites) coupled via electronic 
couplings and weakly interacting with their environment. Each molecule is modelled 
as a two- level system in its electronic coordinate, where the two levels correspond to 
the ground and the first excited electronic states. The electronic degrees of freedom 
are assumed to be diagonally coupled to a bosonic bath of harmonic oscillators. We 
assume that there is one excitation in the system that can be coherently shared between 
multiple molecules (a molecular exciton). The Hamiltonian of this system is given by 

ij i=l k k 

where is the electronic coupling between the ground state of site i and excited state 
of site j, the ket \i) represents the excited state of site i, hik is the electron-phonon 
coupling between site i and mode k of the bath, {b* k , are the raising and lowering 
operators for mode k of the bath, and LOf. are the bath frequencies. 

Due to the assumption of weak system-bath interaction, the exciton states are 
delocalized and the influence of the environment is to induce jumps between them at 
stochastic intervals. In the steady-state the average population of the exciton states is 
given by the Boltzmann factor: Pj = e~P uli /Z, where P, is the population of the exciton 
state with energy cuj and Z is the partition function. From a microscopic standpoint, 
dynamic transfer of excitation still persists in the steady-state. For instance for a system 
of two sites, steady-state implies that fc 1 _ 5 . 2 Pi = k 2 ^P 2 , where k^j is the transfer rate 
from the eigenstate i to the eigenstate j. Moreover, the ratio of the transfer rates satisfies 
detailed balance: ki^, 2 = k 2 ^ie~ /3uJ12 , where u i2 is the energy difference between the two 
states. k 1 ^ 2 P 1 is also the average number of jumps per second between the two states: 
upward and downward jumps occur with equal frequency in the steady-state. If the 
environment is perturbed out of equilibrium, detailed balance no longer holds and the 
frequency of jumps deviates from the thermally equilibrated rate. 

To compute the statistics of jumps in the steady-state, one can compute the 
cumulant generating function of this process directly from the Lindblad equation. The 
Lindblad equation as applied to EET, for weak system-bath coupling with secular 
approximation is given by [T2] 



a(t) = -^iJ,, + ^ 7 ( W ) (i) Al(^) " ^L4( W ),^(t)}] (9) 

u),m 

where the commutator describes the unitary evolution and the second term describes 
relaxation and dephasing in the eigenstates {\a)} of the unperturbed system 
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Hamiltonian H s (neglecting the Lamb-shift Hamiltonian) . The Lindblad operators are 
given by A m {uj) = X)w=e -e , c m( a ) c m( a ')\ a ) ( a '\ where {u} is the energy difference 
associated with single excitation eigenstates, and c* m (a)c rn (ot') is the exciton interference 
at site m. The factor 7(0;) determines the rate of energy transfer between two 
excitons and is related to the spectral density of the bath J(oj), via the expression 
j(co) = 27rJ(|o;|)|n(a;)|, where n(u) = [e^ — l) -1 is the thermal occupation number. 
In the present treatment it is assumed that each molecule interacts with an identical 
and independent phonon bath. The dissipative term in the Lindblad equation describes 
the incoherent jumps between the exciton states as induced by the phonon bath. In 
the quantum jump approach the Hamiltonian evolution of the quantum system is 
interrupted by jumps at stochastic intervals. A sequence of stochastic jumps constitutes 
a trajectory of the system, and an averaging over a large number of trajectories recovers 
the reduced density matrix a(t) [22] . 

If K jumps are observed after a time t, the (unnormalized) reduced density matrix 
of the system prior to detection is given by the projection of the total density matrix onto 
the subspace of K events [H [10]. We denote this projected density matrix by a^ K \t). 
Conversely, the total density matrix can be written as an expansion of all projected 
density matrices, that is a(t) = YlK=o a The probability of observing K events 

after a time t is given by Pt{K) =Tr[cr^(t)]. At this point for the sake of clarity we 
restrict the analysis to the statistics of jumps between two exciton states separated in 
energy by Note that the total system may still contain an arbitrary number of 
excitons. It can be shown that the Laplace transform of the projected density operator, 
defined by the expression a s (t) = a^e~ sK obeys the equation d s (t) = W s [a s ], 
where the superoperator W s is defined to be 

W s [a s (t)} = -i[H s ,a s (t)\ +^2-y{u} 1 )e-'A m {u} 1 )a B {t)Al{u 1 ) (10) 

m 

+ J2 [ 7H4(w)(7,(i)4,(w) - -J]7(w){4»(w)4nM,<7.(t)} • 

Note that a s (t) has the same structure as the variable qt(K,s) in Eq. ([7]). Typical 
(equilibrium) trajectories of a s (t) thus correspond to rare trajectories of the original 
master equation. The ensemble of trajectories generated by Eq. (ITT]) is known as the 
s-ensemble [13] . Physical time evolution is retrieved when the activation parameter is 
zero (s = 0), whereas s^O contains information with regards to rare trajectories of the 
system, s < (s > ) generates the active (inactive) trajectories where the number of 
jumps are significantly higher (lower) than average. 

The Markovian nature of the Lindblad equation enables us to compute the 
statistical moments of the number of jumps in the steady-state via the large deviation 
principle. Large deviation theory states that for a Markov process, the largest real 
eigenvalue of the superoperator W s is the cumulant generating function of the statistics 
[2]. Following the notation of the previous section, we denote this generator by 0(s). 
In the next section we compute the cumulant generating function and the associated 
statistical moments for a system of two, three, four and seven sites. 
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J(») = ^^- 2 (11) 



4. Results and discussion 

Adopting the electronic and vibrational parameters of the FMO complex, we first 
consider two and three sites (chromophores) and study the statistics of trajectories that 
contain jumps between two given exciton states. We subsequently include further sites 
in the system and study the statistics of selected jumps. The coupling of the electronic 
system to the vibrational environment is described via the spectral density 

2E r ujoj c 
IX UJ 2 + uj' z c 

where E r = 35 cm -1 is the reorganization energy E r = J °° ^p-du, and lu c = 150 cm -1 
is the cutoff frequency of the bath [231 121] • We compute the matrix W s and solve for 
the cumulant generating function 6{s). Figure [T] (a) and 1(b) illustrate the exciton 
states where two and three sites are included in the dynamics respectively. The two site 
model includes sites 1 and 2, and the three site model includes sites 1, 2 and 3. The 
site indices are the standard indices assigned to the chromophores in the FMO complex 
[25] . The exciton states are labelled via a greek letter and a number, where the greek 
letter represents the position of the exciton state in the energy ladder, and the number 
indicates the site that carries the maximum exciton amplitude. For instance \a : 1) 
represents the lowest eigenstate, and indicates that this eigenstate is predominately 
localized on site 1. Due to significant exciton localization in the FMO protein, each 
exciton can be unambiguously assigned to a unique site index. 

Figure [1] (c) and (d) illustrate the cumulant generating function versus the 
activation parameter s for different temperatures. The plots suggest that the two cases 
have similar statistics in the active phase, but begin to diverge close to s — 0. Figured] 
(e) illustrates the Mandel parameter for the case of two sites. Mandel parameter is an 
indicator of the nature of first-order statistical correlations in the system and is defined 
to be 

(k) 0'(s) 

Q = implies Poissonian statistics where probabilities of subsequent events are 
uncorrelated. Q > (Q < 0) indicates super-Poissonian (sub-Poissonian) statistics 
where subsequent jumps exhibit first-order correlations. 

For the case of two sites, the Mandel parameter is negative for all activations. 
This indicates sub-Poissonian correlations (anti-bunching) where a waiting period for 
re-excitation is required before the system can relax back to the lower state. As the 
activation parameter is increased, the Mandel parameter approaches zero (Poissonian 
limit). In this limit there are so few jumps that correlations vanish asymptotically. 



It is striking that anti-bunching is observed even though the transitions are 
thermally activated, and the steady-state dynamics obey a classical rate equation. The 
classical nature of the correlations can be confirmed by considering the classical rate 
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equation characterized by the following dynamical matrix 



W dass (s) 



- H 
K 



re~ s 

-r 



where k and T are the relevant equilibrium transfer rates, and V = ne 1 as required 
by detailed balance. By diagonalizing W c i ass (s), we arrive at the following expression 
for Q(s) 

-2kV 



Q(s) 



< Vs 



(13) 



(k + T) 2 -AKT(l-e- s ) 

Indeed, the same expression for Q(s) is obtained if the full Lindblad equation with de- 
phasing terms is instead considered. The formalism thus captures the fluctuations of 
a quantum system due to its coupling to a non- equilibrium bath. For the system under 
consideration the nature of these fluctuations is entirely classical. 



The corresponding plot of the Mandel parameter for a system of three sites is shown 
in Fig. [1] (f). The trajectories exhibit a dynamical crossover in their statistics, i.e.: the 
correlations change from sub-Poissonian (Q(s) < 0) to super-Poissonian (Q(s) > 0) 
as the activation parameter is varied. The crossover point is temperature dependent 
and is closer to s = at higher temperatures. Transitions close to s = are more 
significant as they indicate that small perturbations from equilibrium can significantly 
change the nature of statistical correlation. Typical trajectories (s = 0) are bunched, 
and the Mandel parameter has a peak close to s — 0. 

One can continue to add more sites to the model and investigate the statistics of 
a particular jump, or the collective statistics of a number of jumps. For instance the 
collective statistics of all jumps to the lowest eigenstate can be investigated. In general, 
beyond a few sites there are many possible jumps and there is no obvious advantageous 
collective counting strategy. The statistics of selected transitions for the four site and 
the seven site models, are explored in the remaining part of this section. 

Four interacting sites exhibit the combination of the features of the two site and 
the three site models, as illustrated in Fig. |2j Figure [2] (b) exhibits the new feature 
that typical trajectories (s = 0) have different statistics at different temperatures: anti- 
bunching at higher temperatures and bunching at lower temperatures. In Fig. [2] (c) the 
trajectories remain anti-bunched over the computed temperature range, similar to the 
results of the two site model. 

We consider all possible pairwise jumps of the four site model and conclude that 
for sufficiently high activations all trajectories become anti-bunched and approach a 
temperature independent but transition dependent value of the Mandel parameter. For 
sufficiently low activations, on the other hand, all trajectories approach the Poissonian 
limit. Trajectories that undergo a statistical crossover, approach the Poissonian limit 
from the positive values of Q{s); those that do not, remain in the Q(s) < portion. 
Approach to the Poissonian limit from the Q(s) > side, is therefore a sufficient 
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Figure 1. a) Excition states with only two sites present and the corresponding 
schematics of the interacting sites. The circle containing both sites indicates that 
they are strongly interacting, b) Excition states in the three site model. The three 
site model consists of an strongly interacting pair, weakly interacting with a third site. 
For each case the statistics of jumps between the excitons associated with sites 1 and 
2 are computed, c) 0(s) for the two site model, d) and the three site model, e) Q(s) 
for the two site, f) and the three site model. The site energies are ei = 200 cm -1 , 
e 2 = 320 cm -1 , e 3 = cm -1 and the electronic couplings are J\i = —87.7 cm -1 , 
Ji 3 = 5.5 cm" 1 , J 23 = 30.8 cm" 1 [24]. 



indicator of the existence of statistical crossovers. These conclusions are also observed 
to be valid for larger systems. 

Finally, statistics of selected jumps from the seven site system of the FMO complex 
are illustrated in Fig. [3j We observe that fast transitions are more likely to show 
statistical crossovers. The transfer rate between two exciton states separated in energy 
by to is given by = 7(w) ^ m |c m (a)| 2 |c m (a')| 2 . In this system, the relative 

magnitude of the transfer rates tends to be dominated by the interference contribution 
J2 m l c m(tt)| 2 |c m (a')| 2 ; as the variation of j(co) over the energy range of interest is small. 
In other words, dynamical crossovers are observed in trajectories accounting for jumps 
between strongly interfering eigenstates. The dynamics of the rare trajectories are 
therefore sensitively dependent on the distribution of the electronic couplings and the 
degree of derealization and interference. For exciton states that are close in energy, 
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Figure 2. a) Exciton states for the four site model, and the corresponding schematics 
of the interacting sites. The four site model consist of two groups of strongly interacting 
dimers. b) and c) Mandel parameter versus the activation parameter for selected 
downward jumps in the four site model. The additional site energy is 64 = 110 cm , 
and the additional electronic couplings are J14 = —5.9 cm -1 , J24 = 8.2 cm , 
J34 = -53.5 cm" 1 [24]. 



j(oo) is large, resulting in fast transfer. Statistical crossovers are therefore also observed 
in this regime, despite the small exciton interference. Example of this behaviour can be 
seen in Fig. |3] (a). 

There are important differences between the statistics of quantum jumps in EET, 
and those of the familiar problems of electron and photon counting that merit further 
elaboration. Electron counting is fundamentally different from photodetection in that 
electrons are not destroyed, but are simply counted [I]. Quantum jumps in EET share 
certain characteristics of each example: they are photon-like as a jump corresponds to 
transitions between eigenstates; and electron-like as the detection scheme is assumed 
to be a simple counter that does not interference with the state of the system. Unlike 
photon counting, transitions between eigenstates are thermally activated (they are not 
accompanied by emission of photons), and unlike electron counting, the jumps do not 
generate a current. The insights of photon statistics and current fluctuations should 
therefore only be applied with caution. 

We conclude by noting that the importance of the Markovian assumption in 
the present treatment poses a difficulty for many systems exhibiting EET whose 
environmental interactions are strongly non-Markovian. The formalism can, however, 
be applied to model certain types of non-Markovianity by extending the Hilbert space 
of the system Hamiltonian to contain the environmental modes that capture the strong 
non-Markovian features [261 127]. 



11 



3 
2 

a 1 



-1 



0.2 

-0.2 
-0.4 
-0.5 





(a) 


e:5)->. 6:7) 



a 




a 1 





(b) 




| £:2>->) y:l> 





368 

260 
245 

155 

82 

-30 



t:2) 

e:5) 
6:7) 

Y :l> 

P:4> 

a:3) 



(f) 



Figure 3. a) to d) Mandcl parameter for selected downward transitions. Colours 
indicate the temperature of the bath: red (300 K), green (200 K) and blue (100 K). e) 
Exciton states of the 7 site system and the corresponding schematics of the interacting 
sites. The seven site system consists of three strongly interacting dimcrs, coupled 
weakly to a seventh site. 



5. Conclusion 

Markovian master equations are an established method of modelling the dynamics of 
many systems exhibiting transport properties. These equations are primarily used to 
compute the ensemble averages, but in fact contain more information pertaining to 
higher order statistical moments that cannot be gleaned unless an unravelling formalism 
is first implemented. By introducing an artificial coupling field, or an activation 
parameter, the s-ensemble theory allows the unravelling formalism to be extended to 
quantum trajectories that have been perturbed out of equilibrium. Such an unravelling 
formalism has been applied in the present article to analyze the statistics of Markovian 
excitation transport. 

The non-equilibrium steady-state of a system undergoing EET exhibits surprisingly 
rich statistics: The features that are often associated with quantum correlations such 
as bunching and anti-bunching can be observed in the steady-state statistics of the 
thermally-activated quantum jump trajectories. System size also plays a key role in the 
non-equilibrium statistics as systems of two and three sites reveal qualitatively different 
correlations: while in the former the jumps are always sub-Poissonian, in the latter they 
can undergo statistical crossovers. Transitions in larger systems show a combination 
of these two characteristic behaviours. Moreover, statistical crossovers are found to 
be activated by the electronic coherences: trajectories that account for jumps between 
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two strongly interfering eigenstates are most likely to exhibit statistical crossovers. The 
structure of the spectral density can also influence the statistics by promoting electronic 
transitions of particular frequencies [28] . 

We leave the matter of the potential implications of the observed transitions to 
excitation transport as an intriguing open question. As the quantum jumps are not 
directly observable, their statistics cannot be directly probed. However, if the statistics 
can be mapped to an observable of the system, their dependence on the structure of 
the bath suggests that the observable can be used as a mean to probe the system-bath 
interactions and non-equilibrium perturbations of the environment. 

6. Acknowledgment 

J. P. Garrahan is gratefully acknowledged for numerous stimulating discussions. H. 
H. also acknowledges the hospitality of University College London. The Engineering 
and Physical Research Council (EPSRC) of the UK and the Natural Sciences and 
Engineering Research Council of Canada (NSERC) are also thanked for financial 
support. 

References 

[1] L. S. Levitov and G. B. Lesovik. JETP Lett., 58:230, 1993. 

[2] L. S. Levitov, H. W. Lee, and G. B. Lesovik. J. Maths. Phys., 37:10, 1996. 

[3] Ya. M. Blanter and M. Buettiker. Shot noise in mesoscopic conductors. Phys. Rep., 336:1, 2000. 
[4] Y. V. Nazarov. Quantum Noise in Mesoscopic physics. Kluwer academic publishers, Dordrecht, 
The Netherlands, 2002. 

[5] D. A. Bagrets and Yu. V. Nazarov . Full counting statistics of charge transfer in coulomb blockade 

systems. Phys. Rev. B, 67:085316, 2003. 
[6] W. Belzig . Full counting statistics of super-poissonian shot noise in multilevel quantum dots. 

Phys. Rev. B, 71:161301, 2005. 
[7] A. Braggio, J. Konig, and R. Fazio . Full counting statistics in strongly interacting systems: 

Non-markovian effects. Phys. Rev. Lett, 96:026805, 2006. 
[8] S. Gustavsson, R. Leturcq, B. Simovic, R. Schleser, T. Ihn, P. Studerus, and K. Ensslin. Counting 

statistics of single electron transport in a quantum dot. Phys. Rev. Lett., 96:076605, 2006. 
[9] Y. J. Zheng and F. H. L. Brown. Single-molecule photon counting statistics via generalized optical 

bloch equations. Phys. Rev. Lett, 90:238305, 2003. 
[10] Y. J. Zheng and F. H. L. Brown. Photon emission from driven single molecules. J. Chem. Phys., 

119:11814, 2003. 

[11] C. Flindt, T. Novotny, and A. -P. Jauho. Full counting statistics of nano-electromechanical systems. 

Europhys. Lett., 69:475, 2005. 
[12] H.-P. Breucr and F. Pertruccione. The theory of open quantum systems. Oxford University Press, 

New York, 2002. 

[13] J. P. Garrahan and I. Lesanovsky . Thermodynamics of quantum jump trajectories. Phys. Rev. 
Lett, 104:160601, 2010. 

[14] H. Touchctte. The large deviation approach to statistical mechanics. Physics Reports, 478:1, 2009. 
[15] A. Budini. Thermodynamics of quantum jump trajectories in systems driven by classical 
fluctuations. Phys. Rev. E, 82:061106, 2010. 



13 



[16] J. P. Garrahan, A. D. Armour, and I. Lesanovsky. Quantum trajectory phase transitions in the 

micromascr. Phys. Rev. E, 84:021115, 2011. 
[17] S. Genway, J. P. Garrahan, I. Lesanovsky, and A. D. Armour. Phase transitions in trajectories of 

a superconducting single-electron transistor coupled to a resonator. Phys. Rev. E, 85:051122, 

2012. 

[18] J. Li, Y. Liu, J. Ping, S.-S. Li, X.-Q. Li and Y. Yan . Large-deviation analysis for counting 

statistics in mesoscopic transport. Phys. Rev. B, 84:115319, 2011. 
[19] J. M. Hickey, S. Genway, I. Lesanovsky, and J. P. Garrahan. Thermodynamics of Quadrature 

Trajectories in Open Quantum Systems. arXiv:1206.5719. 
[20] M. Esposito, U. Harbola, and S. Mukamcl. Noncquilibrium fluctuations, fluctuation theorems, 

and counting statistics in quantum systems. Rev. Mod. Phys., 81:1665, 2009. 
[21] A. Budini. Large deviations of ergodic counting processes: A statistical mechanics approach. 

Phys. Rev. E, 84:011141, 2011. 
[22] M. B. Plcnio and P. L. Knight. The quantum jump approach to dissipativc dynamics in quantum 

optics. Rev. Mol. Phys., 70:101-144, 1998. 
[23] M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru-Guzik. Environment-assisted quantum walks 

in photosynthetic energy transfer. J. Chem. Phys., 129:174106, 2008. 
[24] J. Adolphs and T. Renger. How proteins trigger excitation energy transfer in the FMO complex 

of green sulfur bacteria. Biophys. J., 91:2778-2797, 2006. 
[25] F. Fassioli, A. Nazir, and A. Olaya-Castro. Distribution of entanglement in light-harvesting 

complexes and their quantum efficiency. New J. Phys., 12:8, 2010. 
[26] L. Mazzola, S. Maniscalco, J. Piilo, K.-A. Suominen, and B. M. Garraway. Pseudomodes as an 

effective description of memory: Non-markovian dynamics of two-state systems in structured 

reservoirs. Phys. Rev. A, 80:012104, 2009. 
[27] H.-P. Breuer. Genuine quantum trajectories for non-markovian processes. Phys. Rev. A, 

70:012106, 2004. 

[28] H. Hossein-Nejad, C. Curutchet, A. Kubica, and G. D. Scholes. Delocalization-enhanced long- 
range energy transfer between cryptophyte algae PE545 antenna proteins. J. Phys. Chem. B., 
115:5243-53, 2011. 



